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Abstract. A model is presented that explains the 'torsional oscillation' pattern 
of deviations in the solar rotation rate as a geostrophic flow. The flow is driven by 
temperature variations near the surface due to the enhanced emission of radiation by 
^SJ ' the small scale magnetic field. The model explains the sign of the flow, its amplitude 

, and the fact that the maxima occur near the boundaries of the main activity belts. 

■ The amplitude of the flow decreases with depth from its maximum at the surface 
but penetrates over much of the depth of the convection zone, in agreement with the 
data from helioseismology. It predicts that the flow is axisymmetric only on average, 
and in reality consists of a superposition of circulations around areas of enhanced 
magnetic activity. It must be accompanied by a meridional flow component, which 

■ declines more rapidly with depth. 
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^-_{ , 1. Introduction 

a^ 

' The so-called torsional oscillation (Howard and LaBonte 1980) was 

. discovered as a small time- and latitude-dependent modulation of the 

"fH \ rotational velocity of the Sun as measured from Doppler shifts. In 

Q^' averages of the azimuthal component of the solar surface velocity field 

Q . (synoptic maps) it appears as a band of increased velocity that drifts to- 

I wards the equator during the sunspot cycle, together with the magnetic 

c/3 ■ activity. Its amplitude is of the order ~ 5 m/s, and is most prominent 

^ ' on the equatorward side of the main activity belt. 

The obvious connection of the flow with the solar cycle, and its 
^ . relatively small amplitude, suggest that it may be a secondary effect, 

\ somehow caused by the magnetic fields that are the main manifes- 

tation of the cycle^. One possibility would be that it is driven by 
uncompensated Lorentz forces; this hypothesis has been put forward 
immediately after the discovery of the oscillation (Schiissler, 1981), and 
has been the basis of several subsequent interpretations. An exception is 
Ulrich's (2001), which proposes that a hydrodynamic mode of the solar 
envelope must play an essential role, since the flow is already present at 
activity minimum, when no active regions are present. The energy in 
the oscillation is small, however, compared with the inferred magnetic 



^ If this is the case, the term 'oscillation' would be somewhat misleading, as it 
suggests a cyclic variation due to some restoring force intrinsic to the flow itself. 
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energy of the cycle, hence the suggestion in Uhich et al. (2002) that 
the oscillation is in step with the cycle because it triggers the eruption 

of active regions. 

Instead of being a direct consequence of Lorentz forces, the oscilla- 
tion could be the secondary result of a thermal effect of the cycle's 
magnetic fields. The field wound up by diff'erential rotation in the 
interior of the convection zone could cast a 'thermal shadow' on the 
layers above it by interfering with the efficiency of convection (Spruit, 
1977; Parker, 1987). Oilman (1992) proposes that Coriolis forces acting 
on the downflow created by such a shadow woTild set up a pattern of 
differential rotation. The advantage of this idea is that it produces a 
concentration of the fiow at the boundaries of the active latitude belt, 
as observed. It is not clear if thermal shadows would be of sufficient 
amplitude for this to work, especially if the shadowing takes place 
near the base of the convection zone. In a turbulent diffusion model 
the temperature changes due to shadowing are quite small, because of 
the very high effective thermal conductivity of the convective envelope 
(Spruit, 1977). They are likely to be even smaller (Spruit, 1997) when 
the extreme asymmetry between upward and downward flows due to 
the density stratification of the envelope is taken into account, since 
this causes the upffows to be much more identical in entropy than in 
mixing length models (for results of recent numerical simulations see 
Stein and Nordlund (1998), Asplund et al. (2000). 

The bulk of the flux and energy of the solar magnetic field is be- 
lieved to be located near the base of the convection zone (D'Silva and 
Choudhuri, 1993; Cahgari et al., 1994; Fan et al., 1994). These authors 
have shown that the phenomenology of magnetic fields observed at the 
surface can be quantitatively understood as due to the eruption of the 
loops from this magnetic layer at the base, as in Leighton's (1969) 
classical interpretation of the solar cycle. Models in which the torsional 
oscillation is a consequence of magnetic fields therefore typically also 
place the source of the oscillation in the deeper layers of the convection 
zone. 

With only surface observations available, this prediction could not be 
tested but with the detailed measurements of the Sun's internal rotation 
made possible by helioseismology this has now changed. The torsional 
oscillation has been clearly detected in the variation of the rotation 
profile during the cycle (Woodard and Libbrecht, 1993b; Basu and 
Antia, 1998, Schou et al., 1998). Perhaps surprisingly, the oscillation 
appears to have its largest amplitude at the solar surface. In the most 
detailed results so far (Howe et al., 2000, Vorontsov et al., 2002) the 
pattern can be followed to a depth of about 100 000 km. The systematic 
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decline of the amplitude with depth, however, is hard to reconcile with 
an origin in the deeper layers of the convection zone. 

The new interpretation proposed in this paper takes the apparent 
surface origin seriously. It proposes that the oscillation actually is i) a 
secondary consequence of the cycle's magnetic fields, and ii) is caused 
by the thermal effect which small scale surface magnetic fields of the 
cycle have on their surroundings, not the Lorentz force. As shall be 
demonstrated below, this provides a simple and easily quantifiable 
explanation for the oscillation. It agrees with most of its observed 
properties and makes detailed predictions. 

The key element of the theory is the observation that horizontal 
flows on length scales like those of the oscillation are dominated by the 
Coriolis force, hence approximately geostrophic, like large scale flows in 
the Earth's atmosphere and oceans. 

1.1. Relation of the oscillation to the magnetic field 

The oscillation pattern travels with the main activity belt in a butterfly 
diagram, but its maximum does not coincide with that of the magnetic 
activity. To illustrate the well known relation between the field and the 
flow, I reproduce in Figure 1.1a synoptic map based on the observations 
made at Mt. Wilson from 1986 tiU 1999 (Ulrich 2001). 

Deviations in rotation velocity are clearly associated with the mag- 
netic activity, but avoid the center of the main belt of activity. As the 
lower panel shows, the velocity deviations correlate with the latitudinal 
gradient of the magnetic activity. On the leading (equatorial) side of 
the belts the gradient is strongest, and the flows most pronounced. On 
the trailing side, the flows are of opposite sign. The boundary of the 
activity belt is more diffuse on this side, and the flows are of lower 
amplitude. 

The representation in Fig. 1 is chosen to emphasize the relation 
which I want to draw attention to. It somewhat dccmphasizes a known 
problem with the oscillation, namely the fact that it starts before the 
magnetic field of the new cycle becomes noticeable on synoptic mag- 
netograms. This problem is intrinsic to explanations that assume the 
cycle's magnetic field to be the cause of the zonal flows: the flows are 
quite strong already before the new cycle, as measured by the main 
magnetic activity indicators, sets in. I return to this problem in the 
discussion section. 

The correspondence between velocity and gradient of the (absolute 
value of the) field strength suggests that the flows are geostrophic, 
like the large scale flows in the earth atmosphere and oceans. In these 
flows, there is an approximate balance between the horizontal pres- 
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Figure 1. Top panels: synoptic maps of the absolute value of the magnetic field 
(left) and azimuthal velocity residual (right) measured at Mt. Wilson Observatory 
(Ulrich, 2001). Bottom panel: composite map in which the absolute value of the 
magnetic field is shown in intensity and deviations in rotation velocity in color, with 
blue (red) for faster (slower) rotation than the mean. Color saturation codes for the 
velocity amplitude. Data have been averaged over 3 rotations. 



sure gradient and the Coriolis force, while the contribution of small 
scale ('turbulent') momentum transport processes is small on these 
length scales. The flow speed is therefore proportional to the hori- 
zontal pressure gradient. The direction of the flow is perpendicular to 
the horizontal pressure gradient, leading to the characteristic cyclonic 
and anticyclonic flows circling around low and high pressure systems, 
respectively (for a systematic treatment see Pedlosky, 1982). 

On the Sun, the pressure gradients driving geostrophic flows could 
be thermal, as in the geophysical case, or could have a magnetic origin. 
The magnetic field could cause fiows directly, through the Lorentz force. 
A second possibility is an indirect, thermal effect they have on their 
surroundings; an example is the thermal shadow effect mentioned. An- 
other is the well known surface effect of magnetic fields: by inhibiting 
convection (in spots) or facilitating radiative loss (in the small scale 
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field), magnetic fields reduce or enhance the radiative cooling of the 
solar surface. This cff'ect is entirely passive: the energy content of the 
magnetic structures is not affected, they only act as 'valves' changing 
the radiative energy flux. 

These thermal surface effects, though indirect, are not small. The 
excess emission due to small scale magnetic fields persists for as long as 
the magnetic structures are present on the surface. The total amount 
of excess energy emitted at the surface by a small flux tube during it 
life can therefore be much larger than its total magnetic energy content 
(e.g. Spruit et al. 1991). I show in this paper that this thermal effect 
produces flows of the observed nature and magnitude. 

2. Geostrophic flows caused by active regions 

The behavior of slow horizontal flows in a rotating fluid is governed by 
the Rossby number and, if the effects of viscous friction are important, 
the Ekman number. If viscosity can be ignored, the rate of change of 
the flow (left hand side of the Euler equation) can be compared with 
the strength of the Coriolis force. If the characteristic time scale r of 
the flow is long, the Rossby number Ro = l/(i7r) is small, and there is 
an approximate balance between the pressure gradient and the Coriolis 
force. This is called the geostrophic balance; it holds, for example, for 
the large scale flows associated with weather systems in the Earth's 
atmosphere. The latitudinal component of this balance is 

dyp = —2pVx^s\n.A., (1) 

where A is the latitude and x, y are local cartesian coordinates with x 
counted positive in the azimuthal direction (in the direction of rota- 
tion), and y positive towards the North pole. 

In flows on smaller length scales L, friction becomes important as 
well; its effects are measured by the Ekman number, E = h'/{QL'^). 
Friction has to be taken into account when E k, 1. Small length scales 
are present at the boundaries of a large scale flow. The effects of viscos- 
ity in a rotating flow are thus typically observed in the form of Ekman 
layers at the boundaries of large scale flows. 

To obtain typical values for the Rossby and Ekman numbers char- 
acterizing the torsional oscillation, take as the time scale r the time 
for the pattern to change noticeably due to its drift in latitude. This 
would be a time scale of the order 1 yr, so the Rossby number is of the 
order Ro ^ 0.01. Assuming a kinematic viscosity u ^ 10^^ as a measure 
for the turbulent exchange due to granulation, we can determine the 
length scale Le on which friction becomes important by setting E ^ 1. 
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This yields Le ~ 10 Mm. On scales larger than this, the geostrophic 
balance holds; on scales of this order and smaller, the effects of friction 
have to be included. 

2.1. The active region belt: a low-pressure system 

The sign of the torsional oscillation, with an increased rotation rate 
on the leading and a reduced rate on the trailing side of the activity 
belts, is in the same sense as in a low-pressure system in the Earth's 
atmosphere. An interpretation of the oscillation as a geostrophic flow 
thus requires that the activity belts are zones of slightly lower pressure. 

By vertical hydrostatic balance, a temperature perturbation extend- 
ing from a depth z = to z = d (counted positive into the Sun) causes 
a pressure perturbation Sp &t z = given by 



where H is the pressure scale height, and n the number of pressure 
scale heights over which the perturbation extends. 

The observed sign of the flow thus requires that the subsurface 
temperature is lower in an active region. This may be the reason why 
the possibility of a geostrophic flow has not been considered much in 
explanations of the oscillation: the observed higher brightness of active 
regions would seem to argue for the opposite. 

In fact, the brightness of active regions produces just the right sign 
of the temperature perturbation. To see this if one has to suspend for 
a moment the intuitive notion that the brightness of active regions 
must be due to an enhanced temperature imposed from below. The 
effect of small scale magnetic fields on their immediate environment is 
a local increase of the heat flux escaping through the surface (Spruit 
1977). A patch of solar surface with small scale fields has a slightly 
greater emissivity, a lower limb darkening, and it radiates a bit more 
like a black body than the normal solar surface. At any given horizontal 
surface the temperature in the magnetic structure is lower than average. 
That it nevertheless emits more is because the opacity in the structure 
is reduced, by the effect of the magnetic pressure. One result of this is 
that radiation is emitted from geometrically deeper layers than average, 
where the temperature is higher. Another is that when seen near the 
limb, the walls of the depression are seen. Even if those are not hotter 
than typical photospheric temperatures, they arc seen face on hence 
bright, as opposed to the surroundings, which are seen limb-darkened 
('bright wall effect'). 

The greater emissivity implies an enhanced rate of cooling of the gas 
exposed at the surface. As it cools at the surface the gas is swept into 
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the intergranular lanes, so the coohng effect propagates down below 
the surface together with the intergranular downdrafts. The net result 
of the presence of small scale fields at the surface is thus a slight 
temperature reduction below the surface. 

Apart from this theoretical prediction there are observational indi- 
cations pointing in the same direction, most clearly from the p-mode 
frequency shifts during the solar cycle. A slight increase of the fre- 
quencies is observed with increasing activity. The Z-dependence of the 
effect shows that it is caused mostly by a change in the wave prop- 
agation near the surface (Woodard and Libbrccht 1993a). Goldrcich 
et al. (1991) noted that this is opposite to expectation if active re- 
gions are heated from below. Higher temperatures would increase the 
propagation speeds, which acts in the right direction, but the thermal 
expansion of the envelope has to be included as well. This decreases 
the frequencies because of the increased path length, and dominates 
the frequency shifts because the expansion is linear in T, while the 
sound speed only increases as ^/T. 

In addition, a temperature increase, if it originates at a substantial 
depth, would give a different dependence of the frequency shifts on 
the degree / of the modes. These conclusions have been substantiated 
with the much more accurate measurements of p-modc shifts with MDI 
(Dziembowski, Goode and Schou 2001; Dziembowski and Goode, 2002). 
The helioseismic observations thus rule out heating from below as the 
cause of the enhanced brightness in active regions. 

The filling factor of the small scale magnetic fields is small, and 
so are the corresponding temperature changes. To see if the effect is 
sufficient by order of magnitude, consider a simple purely geostrophic 
estimate, based on eq. (1). If L is the length scale of the perturbation 
in the latitudinal direction, the relative pressure perturbation is related 
to the flow speed v by 

6P/P = 2smA^, (3) 

^s 

where Cg = (p/p)^^"^ is the (isothermal) sound speed. With ~ 5 m/s, 
L of the order of the width of the belt of activity, L ~ 3 10^° cm, and 
Cg ~ 10^ cm/s: 

dP/Pr^ 510-^. (4) 

This small pressure effect can be produced by a temperature change of 
the same order or less (cf eq. 2). 

This estimate of the flow speed does not take into account the ef- 
fect of viscosity, however. The pressure difference needed for a given 
velocity will turn out the be significantly higher if friction is taken into 
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account. This is done in the more detailed model developed below. For 
the moment, the estimate suffices to show that the temperature and 
pressure effects associated with geostrophic flows on the Sun can be 
quite small. 



3. Quantitative model 

To quantify these ideas better, a model is needed for the depth de- 
pendence of the temperature perturbations, as well as a model for the 
stratification of the convection zone. 

3.1. Stratification 

The density in the convection zone varies over several orders of magni- 
tude over the range of depths of interest here (several 10 Mm), hence 
its stratification has to be taken into account for a meaningful quan- 
titative estimate. If it were just the stratification of density that is 
of importance, a simple polytropic stratification would be appropriate. 
But since the flows involve differences in temperature, the stratiflcation 
of entropy has to be treated a bit more realistically than just using a 
simple isentrope. I adopt here the 'pseudopolytropic' model used before 
in Spruit (1982, 1991). As in a standard polytrope, the temperature in 
this approximation varies linearly with depth into the Sun: 

r = roC, c = i + ^Ao, (5) 

where z = at the photosphere, and zq a fitting constant to be de- 
termined. A good fit to the stratification of the mean density in a 
mixing-length model of the convection zone is 

P = PoC', (6) 

if zo is chosen of the order 

zo~ 4.510^. (7) 

The gas pressure 

P = Po^, (8) 

is in hydrostatic equilibrium with Pq = zopog/S. The logarithmic tem- 
perature gradient is 

This stratification corresponds to a polytropic index n = 2, instead of 
the value 1.5 for the isentropic stratification of a fully ionized gas. The 
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value n = 2 gives a better overall fit to the solar convection zone. This 
reflects the fact that ionization is partial over a significant fraction of 
its depth, in particular the first 10-50 Mm that are important here. 

Also because of partial ionization, the adiabatic gradient Va is not a 
constant in the convection zone. The 'pseudo' part of the pseudopoly- 
tropic model now consists of manipulating Va a bit. I make it a function 
of depth, such that the superadiabatic gradient V — Va has a depth 
dependence which approximately fits the mean stratification of entropy 
in the envelope, 

G = V-Va = GoC"^ with Go -0.25. (10) 

In this way, a realistic entropy gradient can be included without giving 
up the algebraic convenience of a polytropic model. 

3.2. The temperature perturbation below the surface 

From the observed emissivity changes at the surface we can estimate 
the enhanced rate of cooling in active regions. In order to calculate 
what this implies in terms of temperature changes below the surface, 
a model is needed for the downward propagation of the cooling effect 
in the intergranular lanes. This depends on the details of mixing and 
entrainment between downdrafts and upflows. These are known in suf- 
ficient detail from numerical simulations of the upper few Mm of the 
convection zone (Stein and Nordlund 1998, Asplund et al. 2000), but 
are increasingly uncertain with increasing depth. 

These same mixing processes, however, also determine the mean 
stratification of entropy below the surface. This mean stratification is 
well determined. In numerical simulations for example, one finds that, 
though the downdrafts generated by cooling at the surface vary greatly 
in their value of the temperature contrast ST{z), the average stratifi- 
cation of entropy is in fact reasonably close to that of a mixing-length 
model (Rosenthal et al. 1999). This stratification is also well determined 
observationally from helioseismology. It now turns out that the mean 
stratification of entropy is all that is needed to compute the effect that a 
small change in surface cooling has on the mean temperatures below the 
surface. This removes the uncertainty associated with mixing between 
upflows and downdrafts. To show this I consider two models for the 
mixing process. In the first, the scenario suggested by the numerical 
simulations mentioned is used. The result is then shown to be almost 
equivalent when a mixing length or turbulent diffusion model is used. 

A downdraft generated as an intergranular lane at the surface with 
an entropy Si is compressed to a smaller volume by the rapidly increas- 
ing pressure. This is offset by mixing with the nearly isentropic upflows 
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surrounding it, which increases the filUng factor of the downdrafts 
again. In the process, however, the entropy difference with the sur- 
roundings decreases. The result is a mean entropy S{z) which increases 
with depth in such a way that the difference between this mean and 
the entropy of the upflows, AS{z) = S{z) — So, decreases with depth 
(producing the so-called 'superadiabatic' layer). 

Horizontal pressure equilibrium is a good approximation in the down- 
drafts, and the entropy difference AS' at a given pressure is approxi- 
mately proportional to AT/T (exactly so for an ideal gas with constant 
7). This difference AT/T between the mean temperature and that of 
the adiabatic upflows is a known function of depth (for example from 
helioseismology). In our pseudopolytropic approximation, it is found 
from (10) by integration: 



In order to find the depth dependence of a small perturbation caused 
by a small change in the mean temperature of the downdrafts, I now 
assume linearity in the perturbation. Thus I neglect the effect of this 
small change on the dynamics of the downdrafts, i.e. treat them as 
passively carried with the downfiows. The depth dependence of the 
temperature perturbation is then given by (11): 



where e = {5T/T)q is the value just below the photosphere (optical 
depth of a few, such that the effects of radiative transfer are already 
small). 

In a conventional mixing length model, the depth dependence of 
the superadiabatic gradient is governed by the mixing length-to scale 
height ratio and the heat flux. By linearizing the dependence of the 
heat flux on the superadiabatic gradient in this model, the effect of 
a small change in heat flux can be calculated. If horizontal turbulent 
transport is neglected, the superadiabatic gradient would just increase 
by the same factor at all depths in response to an increase in heat 
flux. One verifies that this would, to within a numerical factor of order 
unity, yield the same result as (12). The main differerence thus lies in 
a possible contribution of horizontal diffusion of heat. Since the width 
of the active latitude belt is of the same order as the depth of the 
convection zone, the effect of such horizontal diffusion will also add 
only a factor of order unity. 

The depth dependence of the temperature effect due to enhanced 
surface emissivity is thus rather independent of the specific model for 
convective heat transport used. With sufficiently detailed helioseismic 




(11) 



6T/T = eC 



-2 



(12) 
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observations of the flow, however, it might be possible to distinguish 
between the predictions made by difi^erent models. 

3.3. First estimate: Purely geostrophic approximation 

With this determination of the temperature effect, the expected geo- 
strophic flow amplitude can be computed as a function of depth. I do 
this flrst in the inviscid approximation; this is then improved in the 
next section. 

By vertical hydrostatic equilibrium^ we get the pressure perturba- 
tion associated with ST: 

dlnP = J ST/Tdz/H = ^^C^ (13) 

where Hq is the surface pressure scale height Hq ~ 1.5 lO^cm. Prom the 
geostrophic balance (3), and the estimate for e given in 2, this yields 
the azimuthal velocity amplitude of the geostrophic flow: 

where Cgo ~ 10 km/s is the isothermal sound speed near the surface. 

Por a quantitative estimate, a plausible value is needed for the am- 
plitude e of the temperature change below the surface. I do this first 
from a theoretical consideration of small flux elements, and then check 
this estimate against the bolometric changes in the solar flux known 
from the measured irradiance variations. 

A small flux clement ('tube') channeling additional radiation through 
the surface (Spruit 1976, 1977) can have a bolometric contrast of some 
10%. This is the total effect, including the area of reduced temperature 
('dark ring') immediately surrounding it, and normalized by the cross 
section of the element. Very small elements can have significantly larger 
values, while in bigger ones (knots, pores) the contrast is actually neg- 
ative, as in sunspots. The value of 10% agrees with observations of the 
smallest elements where this contrast can be measured directly. An a 
priori estimate of the effect of the small scale field in the activity belt 
thus depends on the distribution of sizes of magnetic elements, which 
is poorly known (see however Spruit and Zwaan, 1981). This is at the 
moment the main uncertainly in comparing the theory to observations. 

^ One may wonder if vertical hydrostatic equilibrium of this temperature per- 
turbation is compatible with the horizontal pressure equilibrium that holds for the 
individual downdrafts. This is a matter of length scale (or Ekman number). On the 
large horizontal scales considered here, the horizontal balance is between pressure 
perturbation and Coriolis force: the geostrophic balance. On the small granulation 
scale, the geostrophic balance breaJjs down, and is replaced by pressure equilibrium. 
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Better known is the filling factor of the small scale magnetic field. 
Averaged over the active latitudes, this is of the order 1%, correspond- 
ing to an average (absolute) surface field strength of 10 G. The 10% 
effective contrast of the flux elements thus translates into a mean 
surface flux increase SF/F ^ 10"^. This can be compared with the 
irradiance variations measured at Earth during the solar cycle (e.g. 
Prohlich and Lean, 1998). These show a peak-to peak variation of about 
7 10~^. Assuming that the main contribution of this signal comes from 
the active latitudes, with a width A A ?s 30°, this corresponds to a 
variation in surface flux within the active latitudes of SF/F ~ 1.5 10~^. 
This is quite comparable to the theoretical estimate above. I assume a 
value SF/F = 10~^ in what follows. 

The number for the difference in cooling rate needs to be translated 
into a temperature change below the surface. The radiation emitted 
by the granulation flow, from the time radiative loss first becomes 
important till the time the flow descends back in the intergranular 
lanes, produces the entropy difference between the downdrafts and the 
upflows. This is known from observations and from numerical simula- 
tions of granulation; it corresponds to a relative temperature change 
(at a given pressure) of order unity (from about 11 000 K to 6000 K). A 
change SF/F in the emitted radiation thus yields a temperature change 
in the downdrafts of the same amount, e ~ SF/F. Thus our estimate 
for the average temperature effect due to flux tubes in the activity belt 
is 

e«10-^ (15) 

Taking for the L the approximate width of the sunspot belt, L 
310^°, and this value for e we get v ^ 5m/s at 2; = 10 Mm. This is 
comparable with the values observed helioseismically at this depth. 

At the surface, however, the estimate predicts much larger values 
than observed, around Vg = 100 m/s. As discussed in 2.1, this is because 
we have neglected friction, which is important over the first 20 Mm or 
so. The surface is dominated by an Ekman layer. This is the subject of 
the next section. 



4. Ekman layer 

In this section we derive the fiow which develops under the influence of 
a steady horizontal pressure perturbation near the surface. The effect of 
mixing by the small scale granulation flows is modeled by an isotropic 
kinematic viscosity coefficient. The resulting model is very similar to 
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the standard model for an Ekman flow in geophysics. Due to the vis- 
cosity, an Ekman layer develops near the surface, which modifies the 
flow compared with a purely geostrophic approximation. 

The main difference with respect to the geophysical case is the 
very strong density stratification in the solar envelope. I take this into 
account using the analytic stratification model presented in section 3.1. 

4.1. Equations 

The equation of motion for steady flow of a viscous fluid in a rotating 
system is 

= -Vp + 2/9V X n + T, (16) 

where r = diver is the viscous force, in terms of the viscous stress tensor 
(jjj. The value of the kinematic viscosity coefficient, from a mixing 
length model for the convection zone, is 

u^hv^ 10^^ cmVs, (17) 

relatively independent of depth. The dynamic viscosity r] = pu \s thus 
a very strongly increasing function of depth. The viscous stress in such 
a stratified fluid is (e.g. Landau and Lifshitz, 1982) 

.dvi dvj. , 2 
'^'^ " ^""^d^j ^ ^ ~ :^i^)Sijdwv, (18) 

where x is the (kinematic) compressional viscosity coefficient. Both ly 
and X are uncertain, depending on the details of the turbulence model 
used, but of similar order of magnitude. Since the essence of the flow 
studied depends mostly on the shear viscosity, I simplify the equations 
by assuming x = so that the compressional viscous terms disappear. 

The effect of the Coriolis forces depends on latitude. For the present 
estimates, it is sufficient to consider the situation near the North pole, 
so that the rotation vector is vertical and pointing out of the Sun. 
Use a local cartesian coordinate system with azimuthal direction x, 
latitudinal direction y and vertical r (counted positive out of the Sun) . 
This vertical coordinate has the opposite sign from that used for the 
stratification model of 3.1, i.e. r = —z. This is done here so that the 
rotation rate $7 will be a positive number. 

The horizontal equations of motion for an axisymmetric flow (inde- 
pendent of x) are then 

— dypP — 2pVx^ + 2pvdyyVy + dr\pu{drVy + OyV.f)] = 0, 

2pVy^ + dr{pi'drVx) = 0, (19) 
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where p is now the relative pressure perturbation 5P/P, and P{r) the 
mean pressure stratification. 

To represent the dependence of the driving pressure gradient on lati- 
tude it is sufficient to take a simple sinusoidal form with wavenumber k. 
The resulting flows then also have a harmonic dependence on latitude 
with this wavenumber, and we can replace all latitudinal gradients by 

dy = ik. (20) 

Further simplification is possible by noting that the depth of the Ekman 
layer, of the order L,. = 10 Mm (sec 2 above), is small compared with 
the latitudinal scale of the flows, of the order ~ 100 Mm or more. 
Expanding eqs. 19 systematically in the small parameter Lr/Ly, we get 

—ikpP— 2pVx^ + dr{pvdrVy) = 0, 

2pVyQ. + dr{pvdrVx) = 0. (21) 

Due to this expansion in L^/Ly, the vertical velocity has disappeared 
from the equations. Hence neither the vertical equation of motion nor 
the continuity equation are needed to solve the problem (if needed, the 
vertical velocity can be reconstructed afterwards from the solution by 
using the continuity equation) . 

Four boundary conditions are needed. Since the driving force van- 
ishes with depth (eq. 13 above), the velocity also vanishes at large 
depths: 

Vy,Vx^Q {z^oo). (22) 

In the atmosphere, the granulation velocity and hence the viscosity van- 
ishes with increasing height. This implies that the appropriate bound- 
ary conditions at the surface 2 = are stress-free: 

drVy = drVx = (^ = 0). (23) 

To see this, approximate the decrease of the viscosity with height as a 
jump at z = 0. Across this jump, the r — y stress must be continuous 
(because there are no surface forces, only volume forces). Since the 
viscosity vanishes above the jump, this stress vanishes. Just below the 
surface, this implies drVy = 0, because the viscosity is finite there. The 
boundary condition to be used is thus stress-free. 

Using the analytic stratification model from the previous section, 
and the depth dependence of the relative pressure perturbation from 
the model of section 3.1, the equations of motion are, in terms of the 
dimensionless depth coordinate C = 1 + z/zq = 1 — r/ zq: 

-ivgC- Cvx + Ed^iCd^vy) = 0, 

(:\ + EdciC%Vx) = 0. (24) 
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In these equations, 

is now a 'surface value' of the Ekman number, i.e. the value based on the 
length scale zq of the stratification near the surface. Prom the model 
stratification of section 3.1 and with a kinematic viscosity v ~ 10-*"^ 
cm s~^ as before, we have E 800. The constant 

ezo kPo 

(with the dimension of a velocity) is the purely geostrophic velocity 
amplitude introduced in section 3.3. It measures the strength of the 
driving force, in terms of the amplitude e of the surface cooling effect. 

4.2. Solution 

By inspection, one finds that the solutions of the homogenous part (h) 
of eqs (24) which satisfy the BC at C — oo are of the form 

{v.,y)n ~ C"' exp[-^^^(-l ± i)], (27) 

while a particular solution (i) is 

Vyi = 0, Va;i = -iVgC~^. (28) 

Replacing the latitudinal dependence e*^^ by its real equivalent, so 
that 

p~ cos (fey), (29) 

and applying the stress-free boundary conditions at the surface C = 1) 
the solution is found to be 

= sin(A;y)C-'[l - cos(C/A)e-^/^], (30) 

Vy = sin(A:y)C"^ sin(C/A)e-^/^, (31) 

where 

A = (2S)V2 ~ 40. (32) 

The direction of the flow rotates with depth, while its amplitude de- 
creases. It is similar to the classical 'Ekman spiral' as is produced, for 
example, by wind stress on the surface layers of the oceans (for a de- 
tailed discussion see Pedlosky 1982, Ch. 4.10). One difference with the 
geophysical case is the additional factor which is a consequence of 
the density stratification. In addition, since the driving force is present 
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Figure 2. Depth dependence of the Ekman flow. SoUd: azimuthal (zonal) component, 
broken: latitudinal component. 

throughout the volume rather than only at the surface, the flow declines 
more slowly with depth than in a classical Ekman layer. 

Since 1/A is a rather small number, the velocity at the surface = 
1) is approximately 

Vyo ~ ""xo ~ y sm{ky). (33) 

The velocity at the surface thus makes an angle of nearly 45° with 
respect to the direction of the driving force, just as in the classical 
Ekman spiral. 

With the value of the surface cooling parameter e as estimated in 
section 3.3, we have Vg 250m/s, and Vyo = Vxo ~ 6m/s. The solution 
with these parameter values is shown in Fig. 2. The oscillatory behavior 

typical of an Ekman layer is still present, but much less noticeable than 
in the geophysical case. This is due to the additional factor in the 
depth dependence. 

5. Interpretation and predictions 

The amplitude Vg of the flow at the surface is now smaller than the 
purely geostrophic estimate (14) by the factor 1/A ^ 0.025. For the 
same value of the driving parameter e it 6 m/s, in excellent 

agreement with the surface Doppler measurements (e.g. Ulrich 2001). 
The model predicts a velocity amplitude that declines smoothly with 
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depth, to about 1/3 of its surface amplitude at a depth 2; 60 Mm. At 
the base of the convection zone, the amphtude is still about 10% of the 
surface value. The flow thus extends over much of the convection zone, 
though with reduced amplitude compared to the surface, in agreement 
with the helioseismic observations (Howe et al., 2000; Vorontsov et al. 
2002). 

The velocity pattern produced by the model is in agreement with 
the observation that the flow is concentrated at the boundaries of the 
main activity belt (see Fig. 1). The sign of the flow agrees, provided we 
accept that the small flux elements in an active region acts as 'leaks' by 
which an enhanced heat flux passes through the surface, as discussed 
in section 2.1. The consequence of this leak is additional cooling below 
the surface in active regions, the opposite of what would be the case if 
they were heated from below. 

The model makes a few testable predictions. 

5.1. PREDICTIONS 

5.1.1. Meridional component of the Ekman flow 

A generic feature of an Ekman spiral is the direction of the flow at 
the surface, which is at an angle of nearly 45° to the direction of the 
driving force. In our application, this means that there must be a flow 
directed from the edges towards the center of the main activity belt. 
Its amplitude must the same as that of the torsional oscillation. This 
feature has, in fact, been observed (Komm, Howard and Harvey 1993, 
Komm 1994). Prom data obtained at Kitt Peak these authors find a 
pattern of meridional flow in just this sense, directed towards the center 
of the active latitude zone at an amplitude of about 5 m/s. 

5.1.2. Depth dependence 

A further prediction is that the direction of the flow changes with depth. 
As Fig. 2 shows, the latitudinal velocity component is comparable to the 
zonal component only in a thin layer near the surface. It declines much 
more rapidly than the zonal component and disappears below about 30 
Mm. This can probably be tested with helioseismic data. Some results 
relevant in this context have already been given by Gonzalez Hernandez 
& Patron (2000). 

5.1.3. Nonaxisymmetric structure 

The model also predicts that the flows induced by active region cooling 
are axisymmetric only on average. Activity patches with scales suffi- 
ciently large compared with the Ekman layer depth ze = zqX ^ 20 
Mm, and sufficiently long lived compared with the rotation period. 
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should each have a cyclonic circulation around it (i.e. directed with 

respect to the rotation axis in the same sense as a low-pressure system 
on Earth). These flows should peak in the outer parts of the patch, 
and be accompanied by a flow towards its center with approximately 
the same surface amplitude. With depth, this flow would become more 
nearly geostrophic . On smaller scales, between 10 and 100 Mm, the flow 
speeds produced by the same field strength would be higher because of 
the stronger horizontal gradients, though the smaller scales would also 
make such features harder to detect. 

Some nonaxisymmetric structure in the flow has, in fact, already 
been detected. Ulrich (1998, 2001) shows the existence of long-lived 
azimuthal orders m ;$ 8. 



5.1.4. Speed of downward propagation 

The model predicts that the flows start at the surface, together with 
their cause: the surface cooling due to the small scale magnetic fields. 
They would propagate down together with the downdrafts, so that it 
would take some time for the flows to develop at greater depths. The 
predicted time delay is small, however. The speed of downdrafts is of 
the order 1 km/s at the surface, decreasing with depth due to mixing 
and entrainment. Since the filling factor of the downflows in strongly 
stratified convection is less than the 50% assumed in mixing length 
models, the downdraft speed is larger at all depths than the mixing 
length value. Assuming as a conservative value for the mixing length 
velocity in the lower half of the convection zone of about lOOm/s, the 
time it takes for the effects to propagate to a depth of 100 Mm is less 
than lO^s, or less than about 2 weeks. 

Compared with the activity cycle, the flows are thus set up essen- 
tially instantaneously, and the prediction is that they closely track the 
distribution of small scale fields at the surface. This can be tested in 
principle, but a precise measurement is complicated by the fact that 
the mean rotation in the convection zone is known only with a finite 
accuracy. Since the mean rotation rate varies mostly with latitude, a 
small change in the assumed depth dependence will introduce latitudi- 
nal kinks into originally vertical structures in the pattern of rotation. 
The results of Vorontsov et al. (2002, their Figure 3) in fact show such 
sideways displacements, of both signs (directions of apparent propa- 
gation both upward and downward). The effect appears to be small 
compared with the latitudinal drift of the oscillation with the cycle, 
however. 
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6. Discussion 

I have presented a geostrophic flow model for the torsional oscillation, 
which identifies it as a consequence of enhanced surface cooling by the 
small scale magnetic field in active regions. Its achievements are: 

— The model is quantitative and makes detailed predictions. 

— It explains the sign of the oscillation and its concentration along 
the boundaries of the active latitude belt. 

— Using observed values for the enhanced emissivity of active regions, 
and a standard value for the turbulent viscosity due to granulation, 
it reproduces the correct amplitude of the flow. 

— It explains why the flow has its maximum amplitude at the surface, 
and reproduces the approximate depth dependence as measured 
with helioseismology. 

The accuracy of the model can still be improved with better input 
data. The most critical data it needs are the bolometric brightness 
contrast of active regions, more generally: all small scale magnetic fields 
including ephemeral regions and the network. The surface observations 
required for this currently have limited accuracy, or they rely on proxy 
indicators like Calcium emission. This will probably improve in the 
future, since these measurements are the same as are needed for an 
improved identification of the cause of solar irradiance variations. 

The main difficulty the model faces is that the correct correlation of 
the oscillation with magnetic activity (as seen in magnetograms) breaks 
down around minimum activity, when the flow is already strong but 
little evidence for magnetic fields of the new cycle is seen. This problem 
is addressed below in 6.3. 

6.1. Active regions as cooling agents 

The model is based on the fact that active regions must be areas of lower 
temperature below the surface. This is perhaps somewhat counter to 
the prevailing intuition. One might feel that, since active regions are 
brighter, they should have higher temperatures. This intuition is also 
somewhat rooted in the idea that the energy flux of the convection 
zone comes from below, so that it is natural to associate higher surface 
temperatures with higher heating from below. 

This intuition is misleading, however. In a steady state, the heating 
from below is equal to the cooling at the surface by radiation into 
space, and it is hard to say which of the two is more fundamental in 
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determining the structure of a convective envelope. In a time dependent 
situation, the two are generahy not equal. Which of the two agents, 
heating from below and cooling at the top, provides the most appro- 
priate physical view when considering changes in the convection zone 
depends on the time scale and location of the phenomenon studied. 
On the time scale of stellar evolution, it is the change in conditions 
at the center that determine the changing luminosity of the star. The 
thermal time scale of the star, ^ l(f yr, is very long compared to the 
phenomena of solar activity, however. Even the thermal time scale of 
the convection zone alone, ~ 10^ yr, is very long in this respect. 

Because of this, changes in surface brightness due to variations in 
conditions in the deeper layers of the convection zone, though not 
impossible, arc much harder to bring about than changes in brightness 
caused by surface effects. The darkness of sunspots due to the suppres- 
sion of convection by their magnetic field is a well-known example of 
such a surface effect. The opposite of this effect happens in small scale 
magnetic fields, which act as leaks in the surface through which an 
increased heat flux escapes from the convection zone (Spruit 1977). 
These effects can be viewed simply as modulations of the effective 
emissivity of the solar surface, i.e. changes in the surface heat flux 
at fixed thermal conditions in the bulk of the convection zone. 

The large differences in time scale and energy required for thermal 
changes in the bulk compared with changes at the surface is due to the 
enormous stratification of the density through the convection zone, by 
a factor of 10^. This has no equivalent in the laboratory- and kitchen 
analogs of convecting fluids on which conventional intuition is based. 

6.2. Flows: a proxy for surface temperature variations 

As shown above (4.2) the velocity amplitude produced by a 0.1% change 
in surface emissivity is of the order 5 m/s, for a nominal horizontal 
length scale of 3 10-'^'^ cm. Velocity amplitudes of this magnitude are 
easily detected with current Doppler and helioseismic measurements. 
The brightness changes themselves are much more difficult to observe 
directly. So far, the sensitivity of the order 0.03% needed for such a de- 
tection has been achieved only with disk-integrated measurements like 
the ACRIM irradiance data. Surface flows like the torsional oscillation 
may therefore be a more sensitive way of detecting surface temperature 
changes associated with the solar cycle. 

6.3. The early start of the flow pattern 

A difficulty with the model presented here is the strength of the ob- 
served oscillation at high latitudes, before the beginning of the spot 
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cycle. In synoptic magnetograms, there is very little evidence of mag- 
netic fields at these latitudes and phase in the cycle (e.g. Ulrich, 2001, 
Ulrich et al. 2002). The model presented here predicts that substantial 
magnetic flux is in fact present there, but in a form which somehow 
escapes detection in synoptic magnetograms. 

This is not a new suggestion. Significant coronal activity at high 
latitudes before the beginning of the spot cycle has been noticed already 
in the 60's, in particular in the Pic du Midi observations of Trellis 
et al. (cf. Leroy and Trellis, 1974). It has also been consistently seen 
in the coronographic observations at Sacramento Peak (Altrock, 1986; 
Guhathakurta et al., 1993). Their possible connection with the torsional 
oscillation was first proposed by Wilson et al. (1987). 

The p-mode frequency shifts during the activity cycle show a similar 
effect: the latitude dependence of the sources derived from these shifts 
agrees well with the surface distribution of magnetic fields (Antia et 
al. 2001), except at high latitude, where the effect is stronger than 
expected from magnetograms (Moreno-Insertis and Solanki, 2000). The 
inferred high-latitude magnetic fields peak around activity minimum 
(Dziembowski and Goode, 2002; Goode et al., 2002). 

The form in which these high latitude fields can exist is signifi- 
cantly constrained by their absence in synoptic magnetograms. They 
would have to exist in the form of small scale fiux with polarities 
nearly balancing within the pixel scale of the synoptic maps. High- 
resolution observations at the relevant latitudes around solar minimum 
could test this possibility. Ephemeral active regions, known to have a 
wider latitude distribution than regular active regions and to be most 
prominent at activity minimum (Harvey et al. 1986, Harvey 1989), 
may be a manifestation of these fields (Wilson et al. 1987). It is also 
possible that magnetograms are less sensitive to magnetic fiux in very 
small scale form, due to radiative transfer effects (e.g. Sanchez Almeida 
2000). 
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